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The problem of electromagnetic field quantization is us uaUy 
considered in text-books under the assumption that the c a 
cupies some empty box . The case when ^ , 4 The auan- 

tization^^ e ^o/^ r the m /ie™ i^the^edium^consisttng^o/ tu/o^unt/orm 

I Th? th cas d e if,e o?% afbftfatfl^ogfn UMfS 

k 9 Jn-J 0 eTeU’eu l e n d ? yW 

SSSrS! CBS* 0/ nqn-uni%rm and tHB-fPM*"* media . Ear U 
TonsidihronlTapTrox^ 

for field operltTs in the case of small poiariratton oftheme- 

d s\Zu °the fafe°fhen d non-un ifirTumf-deUdlni dielectric medium 
is clnfifed in some space region with time-dependent boundaries. 

The basis of the subsequent consideration is \^e system of 
Maxwell's equations in linear passive time-dependent dielectric 
and magnetic medium without sources . 


rot E = — i/c dB/dt, 
div D = 0. 

D = s(r,t)E, 


rot H = 1/c dD/dt, 
div B = 0. 

5 = 


(1) 


Introducing the vector potential according to the relations: 

B = rot A, E= -t'o dA/dt (2) 

and imposing gauge conditions 

div(edA/dt ) = 0, <P = 0 (3) 

we can replace the system of the first-order equations (1) with 
the single second-order one: 

rota/t* rotA ) = -1/c 2 d/dt(eM/dt) = He dD/dt. W 

One can check that the vector equation (A) coincides with the set 
of Euler's equations 
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(5) 


d/dt dL/d(djA t ) + d/dx i dL/did^) - dL/dA b = 0 
for the Lagrangian density 

L = 1/2 [ sir ,t)(dA/dt) 2 /c 2 - (rotA) 2 /p{r ,t) ] (6) 

in the case of quite arbitrary time and space dependences of the 
dielectric and magnetic permittivities. Then introducing the ca- 
nonically conjugated variable 

P = dL/didjA) = e(r,t)/c 2 dA/dt = -1/c D (?) 

one can construct the Hamiltonian density 

H - P dA/dt - L = 1/2 [ Z? 2 /e(r, t) + B 2 /p(r,t) ] (8) 

which leads again to eq.(k). But in the general case the expres- 
sion (8) is by no means the energy of the system due to possible 
time-dependences of the coefficients . This fact complicates the 
quantization procedure. The usual procedure consists in introdu- 
cing the field expansion over mode functions 

D{r,t ) = I q N (t)u M (r), B{r,t ) = Zp N (t)u N (r) . (9) 

Substituting these expansions into the Hamiltonian density (8) 
and integrating it over the space variables one gets usually {due 
to certain orthogonality properties of functions u H and u M ) a sum 
of independent oscillator-like Hamiltonians 

H = J H(r,t)d 3 r = 1/2 Z(p N p N 2 + n H q N 2 ). (10) 

After this the coefficients p H and q H are proclaimed operators 
satisfying canonical commutation relations , so that the fields 
become quantized. But this sketch of standard quantization scheme 
shows distinctly that it can be used only in the case when the 
solutions of Maxwell' s equations can be factorised into the pro- 
ducts of two functions: one dependent only of time, and another 

dependent on space coordinates only. In the general case of non- 
uniform and time-dependent medium such solutions do not exist, 
and the usual scheme of quantization is impossible . This means in 
particular , that we cannot obtain any Hamiltonian and. con- 
sequently , any unitary evolution operator. Therefore the Schro- 
dinger 'picture does not exist in the general case. But the He- 
isenberg description is still possible. It can be introduced as a 
generalisation of the approach used earlier by Moore {Ref .11) for 
the field quantization 'in the empty space region confined with 
moving boundaries . 

First of all we notice the important property of equation 
(A): it admits a time-independent scalar product of any two dif- 
ferent solutions in the following form: 

(( A,B )) = -1/2 i Jd 3 r e(r , t) [AdB m /dt - B m dA/dt]. (11) 
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It is essential that the dielectric permittivity is a real 
function , i.e. the medium is assumed lossless. Besides, the vec- 
tor potential has to turn into zero at the surfaces confining 
the integration domain. The case of moving boundaries (conside- 
red in £e/.ll) is included to general situation automatically . 

Suppose that before some instant of time ( let it be t= 0) both 
the medium and the boundaries were time-independent . Then solu- 
tions of (4) could be factorized : 

A(r,t ) = g(r)exp(-iut) , (12) 

ro£[l/p(r) rotg] - co 2 /c 2 s(r)g = 0. (13) 

The scalar product (11) was proportional to the usual scalar pro- 
duct: 


(( A,B )) = l/2(o> fl + UB)erp[i(u B -tt fl )i](0B ( g fl ), (14) 

(06,0<i) = tt*r e (r)g b *g a . (15) 


But it is known that solutions of eq.(13) form the complete 
orthogonal with respect to scalar product (15) set of vector 
functions . Therefore any real vector field can be decomposed 
over this set of functions: 

A(r,t ) = I[a H g H (r)exp(-i(j) H t) + a H m g H m (r)exp(io u t) ] . (16) 

Comparing (14) and (15) we can see that the basis functions can 
be normalised in such way that they will satisfy the relations 

((i4 N ,^4 H )) = S N h> ((^mi^m )) = 0. (17) 

After the instant when the properties of the medium became time- 
dependent the basis functions change their explicit expressions, 
but the scalar products (17) will not change, and instead of (16) 
we can write the decomposition 

Air, t) = na H A H (r,t) + a H w A H *(r ,t)] . (18) 

Then we proclaim the ( time- independent ) coefficients of this ex- 
pansion operators satisfying bosonic commutation relations and 
thus obtain the quant izea field from a classical one. 

If in some time the medium will become again time-independent 
then the physical states will be described with monochromatic 
mode functions of the type (12), which will not coincide in gene- 
ral with the basis functions of expansion (18). Therefore we will 
have two different expansions of the field operator: expansion 
(18) over the states corresponding to the physical photons in re- 
mote past and expansion like (16; over trie physical states ari— 
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sing in future. Designating the " physical " states with the super- 
script "zero", we can expand each set of basis functions into 
a series with respect to another one : 

- $[ a NM^N <0> + Pmm^h <0) (19) 

The corresponding expansion of "new" creation and annihilation 
operators over the set of "old" ones is as follows: 

Qh<0> = + PnH (20) 

The initial state of the quantized field was determined with res- 
pect to the set of "old" operators ( without the superscript 
zero" ). Then using expansion (20) we can calculate all quantum 
statistical characteristics of the field in the final state. Tak- 
ing into account conditions (17) and the evident properties of 
the scalar product (11) 

(01,0)) = ((M))" = -((B* ,A*)) (21) 

one can express the coefficients of expansions (19) or (20) as 
follows: 

a NN = ( (>4 n ,j4 h ( 0 ) ) ) , Pnh = ((^m • (22) 

The quantization scheme sketched above can be applied to the 
most general situation of an arbitrary space-time nonuniform me- 
dium and moving boundaries . However, the explicit calculation of 
the mode functions and coefficients of the canonical transforma- 
tion (20) can be performed only for rather simple special cases. 
The first of them corresponds to the media with factorized elect- 
ric and magnetic permittivities: 

e(r,t) = e(r)%(t), p(r,t) = p(r)v(£) (23) 

( the boundaries do not move). Then mode functions can be also 
sought for in a factorized form: 

A(r,t ) = g(r)^(t) D(.r,t) = e(.r)g(r)r\(t) . (24) 

Let us demand the function g(r) to satisfy the following equation 

roHi/p rotg ) = k 2 e(r)g, k = const. (25) 

Then eqs.i 2) and (4) result in the following ordinary differenti- 
al equations for time-dependent factors of the vector potential 
and electric displacement : 

H = fr 2 cf;/v(t) £ = -cr\/%(t) (26) 

£qrs. (26) resemble equations of motion of an oscillator with time- 
dependent mass and frequencies . The role of generalized coordi- 
nate is played by the electric displacement time-dependent fac- 
tor, while the vector potential time-dependent factor plays the 



role of generalised momentum. Eqs. (26 ) can be replaced by the 
following second-order differential equation: 

f|'+ r(On + Q 2 (t)r\ = 0, t = v/v, Q 2 = k 2 c 2 /v(t)x(t) . (27) 

We shall consider the field inside a resonator. Then solutions of 
eq.( 25) can be chosen real vector functions satisfying the ortho- 
gonality conditions: 

$d*r s(r)g K (r)g L (r) = Jr 2 6 KL . (28) 

Complex solutions of eq.( 27) can be normalized as follows: 

vU) Cnn w - n"n] = -21. (29) 

77its means that we choose the solution of eq.( 27) in the statio- 
nary case in the form of 

Ho ( £) = (v 0 Q 0 )' ! ,2 exp(-i%t) . (30) 

Due to (28) coefficients (22) are not equal to zero only for 
coinciding indices ( intermode interactions are absent) .so we can 
omit the indices. Taking into account eqs. (11) , (26) , (30) one can 
represent these coefficients as follows: 

a = l/2(v 0 /Q 0 ) 1 / 2 (Q 0 q + ir\)exp(iQ 0 t) (31) 

p = l/2(v 0 /Q 0 ) 1 /2 (£ 0 r| - ir\)exp(iQ 0 t) . (32) 

Let us introduce the quadrature components and their variances as 
follows: 

X i = (2)" 1 / 2 (a 0 +a 0 *) * 2 = i(2)' 1 / 2 (a 0 * - a 0 ) (33) 

o I3 = l/ 2 <J t J 3 +X 3 X l > - <X I ><X 3 >. (34) 

Suppose for simplicity that initially the field was in the cohe- 
rent quantum state. Taking into account eq.( 20) one can easily 
obtain the following expressions: 

0,4 = i/2la+pl 2 = i/2v 0 Q 0 lr|l 2 , (35) 

o 22 = l/2la-pi 2 = l/ 2 v 0 Q 0 _ 1 Iql 2 , (36) 

o 12 = Im(ap") = l/2v 0 /?e(qn“) - (37) 

We see that time-dependent medium transforms an initially cohe- 
rent state to a correlated quantum state characterized by 
a nonzero covariance (37) and unequal variances (35) and (36). 

This state minimizes the generalized uncertainty relation by 
Schrodinger and Robertson (ReTsA2 f i3): 
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> 1/4 


(38) 


° H °22 °12 

( the equality takes place in the case under study due to eq. 
(29)). For the detailed review of various forms of uncertainly 
relations see (Ref .14) . Properties of correlated and squeezed 
quantum states were investigated in (Refs. 15-19) . 

Let us consider as an example the case of a parametric exci- 
tation when the properties of the medium harmonically oscillate 
with twice frequency with respect to some ( resonance ) mode. This 
can be achieved, for example , by means of changing the density of 
the medium. Since the magnetic effects are extremely weak, we can 
write 


Q 2 (t) = Q 0 2 (1 + <ecos2Q 0 t) , r = 0. (39) 

We look for the solution of eq.( 27) in the form 

H (t) = (v 0 Q 0 )" 1/2 [ u(t)exp(iQ 0 t ) + u(t)exp(-iQ 0 t)] (40) 

with slowly varying time dependent amplitudes . Substituting (39) 
and (40) into X27) and performing averaging over fast oscillati- 
ons we arrive at the equations (neglecting the second order deri- 
vatives of slowly varying amplitudes) 

u = iQ 0 aeu/4, v = -iQ 0 «u/4 (41) 

whose solutions are 

u(t ) = cosMQ 0 ae£/4) , v(t) = -i sinh(Q 0 xt/h) . (42) 

The variances (35), (36) oscillate with the twice resonance fre- 
quency, but their ratio (the so called squeezing coefficient ) is 
confined at every instant between the values 

exp(-Q 0 <et) < Q 0 2 o 11 /o 22 < exp(Q 0 &t) . (43) 

Certain inequalities for the squeezing coefficients can be found 
for arbitrary time dependence or the frequency in eq.( 27). Consi- 
dering again nonmagnetic medium one can prove the unequalities 
(Refs 19,20) 

[(l-/? 1 /2 )/(l+/? 1 /2 )] 2 <Q 0 2 o 11 /o 22 < [(i+/? 1/2 )/(l-tf 1/2 )] 2 (44) 

where R is the energy reflection coefficient from the effective 
potential barrier corresponding to eq.( 27). 

We would like to emphasize once again that we have used the 
Heisenberg picture for the description of the quantized electro- 
magnetic field. However, since for a factorizable medium (23) 
Maxwell's equations can be derived from the Hamiltonian (8) or 
(10). the Schrodinger description is also possible in this case. 
We shall illustrate it for the most simple case when 


204 


I 



(45) 


v(0 = 1, % = %(ht) , x(°) - 1- 


Here h is a characteristic frequency of the medium properties 
changing. The method used below can be easily applied to a more 
general case. Let c = i and the dimensionless time t 0 = kt. So 
127) results in 


r) + o) 2 ($t 0 )n = 0 , 


(46) 


where 


u 2 (a it 0 ) = Q 2 (ht) = i/x(ht) , 


* = h/k. 


The quantization of a harmonic oscillator with a variable 
^ uency is done by introducing integrals of motion operators 

a(t 0 ) = i(2)' t/2 (y(t 0 >£ - y(t 0 )r\), a*(t 0 > = CaCt 0 )3*. 

//ere y(i 0 ) a " ruling solution " o/ eq.(46) satisfying the 
independent condition 

y(i 0 )y*U 0 ) - y"(£ 0 )y(i 0 ) = 2i 

to ensure the following commutation relation 


(47) 

fre- 

(Ref. 

(48) 
£ime- 

(49) 


[a(£ 0 ) , a + ( t o ) ] = 1 . (50) 

We want to stress the difference in signs of right-hand sides of 
(29) and (49) due to the difference between Heisenberg and bchro- 
dinger pictures. For the integrals of motion to coincide at t 0 -0 
with creation and annihilation operators, it is necessary (in ac- 
cordance with (49)) to choose the initial conditions 


y( 0) = 1, y (0) = i. (51) 

After the instant t Q = t F , when the medium properties stop chang- 
ing, "new" creation and annihilation operators are to be introdu- 
ced: 


a (0) = i(2o )(t F ))' ,/2 (£ - iq «(t F )). (52) 

Then the expansion of "new" operators over "old" ones is 

a (0) = a a(t F ) + pa*(t F ) (53) 

and we can get 

p(t F ) = 1/2 (co(t F ))- ,/2 (uU F )y(t F ) + iy(£ r )). (54) 

Creation and annihilation operators mixing results in a change of 
occupation numbers in a given mode. Generation of photons from 
vacuum due to the medium properties change is worth considering . 
If at t 0 = 0 the number of photons was zero then at an instant t 0 
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>0 the average number of photons is 
n = p"U 0 )p(t 0 ). 

Yablonovitch (Ref. 21) stated that photons will have a thermal 
distribution at a temperature proportional to the rate of the me- 
dium properties change , i.e. proportional to h. Then for high 
photon energies ( proportional to k) there will be 

h ~ ex p(-const k/h) . (56) 


The solution of (39) at a little t 0 can be found for an arbitrary 
law of change of %(t). At t 0 = 0 we shall have 

d/dt 0 ((u) 1/2 p) = 1/2 cfc>((*t 0 ))/df 0 (57) 

p - t 0 /2 d(o(*t 0 ))/dt 0 (58) 

t -> 0 

n ~ (t/2 d(o(ht))/dt). (59) 


For any instant t a decomposition of the solution can be found in 
the limit ae -> 0 ( for high photon frequencies ) bu means of the 
method of multu-scale asymptotic decompositions (Ref .22) for an 
arbitrary law of change: 

y = (<o(a ei 0 ))" 1/2 { exp(is) + ae [ E(&t 0 )exp(is ) + 

+ F(dit 0 )exp(-is)] + ... } (60) 

where 

s = 2/ae [exp($t 0 / 2) - 1] (61) 

and functions E and F are determined by 
F = const 

, , , (62) 
dE/d(st 0 ) = -i/(4« 2 ) d 2 <j)/d(&t 0 ) 2 + 

+ t3/(8<o 3 ) [d<j)/d(&t 0 )] 2 

with initial conditions 

E( 0) + F( 0) = 0, i (E( 0) - F( 0)) = 1/2 d<o/d(aet 0 ). (63) 

For the case ^ = exp(-ht ) , which approximately describes dielect- 
ric permittivity falling achievable in experiments (Ref. 21), one 
can find the exact solution - a linear combination of two Han- 

kel's functions H 0 (i) [ 2/ae exp(ae£ 0 /2)3 and H 0 <2) [ 2/ae exp($t 0 /2) ] . 
Nevertheless an asymptotic decomposition is still more useful and 
can be expressed in an explicit form: 

y = exp(-&t 0 /k)exp(is) + ae exp(-&t 0 /A) { (-i/16) x 

x [1 + exp(-*t 0 /2)]exp(is) + i/8 exp(-is) } + ... (65) 
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Now we have 

p = i*/8 [exp(-is) - exp(-xt 0 /2)exp(is )] (66) 

and the number of photons 

h = * 2 /64 [1 + exp(-xt 0 ) - 2cos(2s)exp(-xt 0 /2) ] (67) 

This number oscillates with a growing frequency and a decreasing 
amplitude and in the limit t 0 -> » lends to a constant 

h -> h 2 /6Ak 2 . 

As we can see it is not in agreement with the statement (38) from 
(Ref. 21). The energy of photons in the mode is 

hwk -> h 2 /§hk exp(ht/2 ), t -> «> (69) 

and it grows without a limit. Also at any time the sum of ener- 
gies of all modes diverges. This can be explained bu the fact 
that it is impossible to decrease dielectric permittivity to zero 
for nondispersive media ; thus the assumption that it does not de- 
pend on a frequency is not valid for high frequencies . In the 
other limit x = h/k -> «> we introduce another dimensionless time 
t, = xt Q . The solution expansion over 0 = i/x ( not valid for 
large t { ) is 

y = 1 + Bit x /2 + 0 2 /4 C 1 1 + i — + ... (70) 

and the first term for the number of photons does not depend on 
the mode frequency: 

P = expit t / 2) - i + Bi/2 [ t i exp(t i /2 ) + i - expt i ] (71) 

h = 1/2 [coshiht/ 2) - 1]. (72) 

Another example of time dependent resonator which can be sol- 
ved is an empty resonator with a moving ideal wall. Moore (Ref. 
11) proposea the following complete orthonormal set of solutions 
(in the special case of a single space dimension, i.e. confining 
with the modes with linear polarisation parallel to the wall sur- 
face) : 

A H (x,t) = (4jiAO -1 /2 {exp[-inNR(t-x)] - exp[-iTiNR(t+x)]) (73) 
(c - 1) . 

These functions depend on the solution of the functional equation 
(L(t) is the position of the moving wall , another wall is assumed 
to be at rest) 

R(t+L(t)) = R(t-L(t )) +2. (74) 
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An approximate solution of this equation in the case of the small 
velocities of the wall was found bu Moore (Ref.ii) and later used 
in (Refs. 23,24). However, that solution is not valid in the case 
of parametric resonance , when 

L(t) = L 0 [i + a sin(2Q 0 t)], I a I«1 , Q 0 = n/L 0 (75) 

(the resonance at the lowest resonator eigenfrequency) . The cor- 
responding solution for small values of the percentaqe modulation 
was founcT in (Ref. 25): 


R(V = t,/L 0 a - 

+ i/4Q 0 £ s£n(2Q 0 


a sirK2Q 0 E) 


+ a 2 [sin 2 (Q 0 O + 


(76) 


Eqs . (11) , (22) , (45) result in the following expressions for the 
transformation coefficients (21): 


NM 1= i/2 (M/N) i/2 \ ° expi in[-NR(L 0 x)±Mx] }dx. (77) 
Pn m J c/Xq-1 

The calculations are rather simple for not very large values of 
time, when the second-order correction in (76) remains small. 
Then the following simple formula for variances can be found 
(Ref. 25): 


o 

o 



= l/2exp(±l/2 naN ) , N» 1, 


la/VI«i 


(78) 


where N is the number of semiperiods of wall's vibrations . One 
can check that the maximum squeezing coefficients given by eqs. 
(43 ) and (78) coincide for equal values of percentage modulati- 
on in two different methods of exciting the field via the para- 
metric resonance . 


Now let us consider the long-time asymptotics for the R-func- 
n under the condition e/L n « 1, e r » L n j’ 


(79) 


t ion under the condition e/L 0 « 1, et » L 0 for an arbitrary 
periodic motion of the wall L(t) -= L 0 + £ fit). We shall choose 
the solution of eq.( 74) in the following form\ 

oo 

R(t) = |_ 0 e N ^ N (O. 

Substituting expansion (79) to eq.( 74) we attain 

PO oo 

^_ 0 £ N /? N [t+L 0 +€/ ( O] = 'j__^£ H R il (t-L <i -£f(t)] + 2. 

Developing both sides of this equation into power series we have 


J_ 0 £ ( * ' ( t+L o )e ( t)/K ' 

x (-{)*'£*' ft (t)/K\ + 2. 


CO oo 



It is convenient to use another summation index M=N+K : 



00 

L 




u, (t+i.o) - (-l'>'R,-L"Ut-L a )]fHt)/K'. = 2 


^“0a“"0 

From this equation we obtain the following system of equations 
for the functions R*(t) (M= 0,1,...): 

(80) 


r=o 


[R H . K iK> (t+L< i )-(-l) t R H - ti { * ) (t-L 0 )]f*(t)/K'. = 25 M0 . 


(81) 


Further we consider the simplest law of motion (75) with e - aL 0 
and the frequency o)n = QQ 0 * " e use the Fouriei transformation me 
thod to solve eq.(8 0) : 

7(o>) = $exp(iut)f(t)dt, fit ) = Soxp(-iG>07(o>)dG>/27T f 
1/2jt ^exp(iwx)dt = 5(o>), 

$expiiot)f(t+L 0 )dt = exp(-ioiZ, 0 )7(o>) , 

$exp(io)t )/ ( N * ( Odt = (-i(o) n 7(g)) . 

Then we jet from (80) the following integral equation: 

1/A/! 5du(-iu) N £„_ N (u) [exp(-tuL 0 )-(-l) N exp(iuL 0 ) ] x 
x $(,)-u(sin M (o)Qt)) = 2 ji6(o))6 mo (^2) 

$ Q (sin H Q 0 t ) = $exp(io)t )sin N (o) 0 t)dt = 

= l/(2jt)"‘* I C,’ (-O 3 5 (u+W-2J)u 0 ). (83) 

Substituting expression (83) to eq.(82) u/e con easily make inte- 
gration over u and arrive at the equation 

-7? M (o))2istn(o)L 0 ) = 26 H0 S(o)) - ^ x 

x [o)+(A/-2J)o) 0 ] n ^ m - n [oh-(A/-2J)o) 0 ] x 
x [exp(-i(j)L 0 )-(-l ) N exp( io>I 0 ) ] (~D 0M * 

Taking into account the formula for derivatives of 6-function 
6 <N) (x) = -A/5 <n_1> (x)/x, 
one can easily find the expression for R 0 io) 

K 0 (u) = 2n6'(o>)/tL 0 . (85) 

Then making inverse Four ier-trans format ion we have 

R 0 it ) = t/L 0 . (86) 

To find the long-time asymptotics we will seek for the solution 
of eq.iW in the form of a sum over 6 -functions 
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K m (co) = £ a K M 6 <M> (u+(/f-2f)u 0 ) (87) 

This choice corresponds to the representation of the function 
R(t) m the form of a power series with respect to the parameter 

(et/L 2 ) and neglecting terms like s 3 t M with J>N. Then only terms 
corresponding to AM are significant in eqASh): 

-7? N (co)2is£n(<*>I 0 ) = 25 N0 6(co) + (-l)°cos(wL 0 , x 

x C (omo 0 )R „ . < (o+u Q ) - (a>-<Og)7? N _ 1 (a)-Wp)] . (88) 

Taking into account the expression for /?„(«) (85) we get 

Ei(<*>) = 2 tt(-1)°/2L 0 cotan(uL 0 ) [ (u+G) g )6' (omo g ) - 
- ((o-<j 0 )6 , (g>-<o 0 )] = 2tt(-1) 0 [5' (w+Wg ) - 6' (a>^>g) ]/2L 0 2 . (89) 

With respect to expansion (87) we have 

a/ = 2tt /iL 0 , ag 1 = -a, 1 = 2ti(-1)°/2L 0 2 . 

Substituting expansion (87) we obtain the following recurrence 
relation for Afe2 

a K " = (-l)°tw g /Z NL 0 [ (^-2^-1 )a c M ' 1 - (N-2K+i )a K _ , N ’ 1 ] . (90) 
Let us introduce the notation a c 2H = ia )C 2H , a c 2M + 1 = a K 2M + I , then 
ffd t H = (-1 ) N a[a K N ” 1 (N-2K-i ) - a K _, N_1 (A/-2AT+1) ] , 

a = (-i)°Wg/2L 0 . (91) 

Making inverse Four ier-trans formation of (87) we get GV>2): 

R H (t) = 1/2ji £ a*" Ut)*expat(N-2K)o Q ) 

Tfien.itis east/ to see that we have got the same expression for 
R\\t) that was given in eq.(76). Now we consider the sum 

<Kt) = E e H R H (t) = E e H Ut) ul Z az H expUto a (N-2K))p c /2n (92) 
N=2 t N-2 K = 0 

for p- 1. Taking into account the evident symmetry condition a c N = 
= -a N . K N we obtain 

$ = 2 E z N (-l) n ? d K H sin(u Q (N-2K)t)p t /2n, z = s t, (93) 
A/=2 K=0 

where n=[N/2 ] anti [ ] is ihe entire part of a number. Taking into 

account that (-l) n = (_i) m <n-i>/2 we introduce the notation $ = 
=-2ImF. with the help of the recurrence relation (91) we get the 
differential equation for the function F 

dF/dz = a[pexp(-io Q t ) - exp{iu 9 t)][zdF/ dz - 2 pdF/dp + 
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(95) 


(96) 


+ (za/jiQ)exp(iu Q t) ] 

with the initial condition F(z=0,p) = 0. Its solution is 

F = -(za/:n:G)exp(ici>gt) + 

where ¥ satisfies the equation 

dV/dz = a[pexp(-io g t) - exp(io 0 t) ] [zW/dz - 
- 2 pdV/dp] + (a/7iQ)exp(iu g t) . 

The particular solution of this equation is as follows : 

¥ 0 = - ( 1 / 2 Jt 0) in[pexp(-io>gt)/(pexp(-ia> 0 t) - 
-exp(icogt))] . 

TTien the generai solution of eg. (95) is the sum of and an ar- 
bitrary solution of the uniform equation (95) with 6-0. The uni 
form equation has the first integral : 

c = exp(-i« Q t/2) (zp 1/2 - i/2a in[ (p 1 /2 exp(-io> 0 t/2)- 
-exp(io) g t/2) )/(p* ' 2 exp(-i(*> Q t/2)+exp( i(»>gt/2) ) ] }. 

Then we get the general solution of eg. (95) in th e f orm of V = 
= + /(?), where fiC) is an arbitrary function of the first in- 

tegral . From the condition ¥(z= 0) = 0 we can determine the form 
of the function f: 

fix) - (l/2irt3) in{ [ i+exp(-2aexp( ia>gt/2)x) ] 2 / 

Mexp(-2aexp( iw 0 t/2)x) }. 

>l/ter some algebraic transformations we find the function t 

^ = -(1/2tiO) tn{ 4 pexp(-i(o g t)exp(- 2 az)/[cos(o)gt/ 2 ) + 

+ 2isin(o g t/2)exp(-2az) ] 2 ). 

7aHna into account the first and the second terms of the expan- 
sion of Rit) in a set of et we obtain the final expression: 

Rit) = t/L 0 - (2/nO) Im Inti + 5 + erp(iu g O(l-?)], (97) 


where a notation E — [ ( — 1 ) ^ *(Og£t/Lo ] is introduced. Now we 

can compute some characteristics of the eiec tromagne t te /ceirf in 

a cavity in the presumption et/L „ s »l, £<<(.<, •. (■?£ us evaLu t °i e Jn 
cumber of photons which will be generated in the resonator with 
moving walls in the Long time limit. For this purpose we must 
evaluate the integral (77; 

£ NH = i/2 ( M/N ) u 2 \exp[-iniiM+N)x+Nfix))] , 
where the function fix) is given by (97): 
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fix) = -2/n 0 Im Zntl+fj+expCtJiQxHl-^)] . 

Let us consider the case when Q=2p is an even number , then f«i. 
D “ e *° ex/L 0 << 1 we can consider fix) as a periodical function 

with the period of oscillations T-2/Q. Then 

d i . . _ . * 


f _ - I " " ' ww-riw < hr I 41 lb I l 

NN ?= iM/N) i/2 I [ exp[-iniN+M)2K/Q] T^ixpi-im. x 
Onh) a=0 ' AM 


x iiN±M)x+Nfix))]dx. 


(98) 


Analyzing the structure of fix) we can approximate it bu three 
linear functions as follows : * 

-(1-60) ix-N+i) , N- 1 < x < AM+i/0-6 

fix) = ii/Q-25) ix-N+[-i /Q) / 8 , AM+1/0-6 < x < AM+1/0+6,’ 

-( 1-60) ix-N+i-2/Q ) , AM +1/0+6 < x < AM+2/0, 

6 = 2 £* /2 /nq. 

Then one can easily calculate integrals in (98) and obtain the 
general expression for the coefficients a NH and 0 NM 

a NH = 2(Af//V) W2 (-l) <N - 1 > (N - M >/[ J r(5W-Af)] sin[ i8tlQ-M)n/Q] x 
x exp[ in(N-M) (0-1 )/0] sin[niN-M) ]/sin[niN-M)/Q] , 

P.» = 2«//V) </ ’(-O < "-‘>'“*">/[n(6JV0 + «)] sin[(6M + rt) I r/0] ( x 9) 
x exp[ iniN+M) (0-1 )/0] sin[ji(AM0 l/sintirOiAO/O] . 

Hereafter we consider only the main resonance of 0=2. After 
S foTthe l coe r fnUenT fp^?^ i0nS We 9&t the followin 9 expression 


I^mh 1 2 = With r 2 ) [1-(-1 ) h cos(2A/6jt) ] [l+(-l) M + M ]/ 

/iM+m) 2 . 


( 100 ) 


r } l f mber of photons in the mode with number M we 
” owing 0 auxi liar y g sum* S “" ° Uer N ' Flrst ' let us evaluate the *ol- 

Siz,x) = ^_cosixN)/[Niz+N) 2 ], ( 10 i) 

where x=26tt«1 and z=M/28»i .Then we have 

00 00 

Si z ,x) = J y ^_cosiNx)exp[-yiz+N) ]/N dy = 

= 1/z -1/2 J yexpi-zy)lni2coshy-2cosx ) dy. (102) 

Since x«l and t/ie region of integration u&i/z«i we can 

? E an * c °/ h yS nd > cosx P° wer series of y anf x Thus up to 

the second order terms we have H 
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S(z,x) = i/z 3 -1/2 J yexp(-zy)ln(y 2 +x 2 ) dy. (103) 

The last integral can be easily evaluated, if one takes into ac- 
count the inequalities y$i/z«x«i. Therefore 

S(z,x ) = -lnx/z 2 + 0(z~ 3 ). (10A) 


The similar sum ( see (100)) 

^ cos(Nx)/[N(z+N)] (-i) M 

can be obtained from (102) bu means of the replacement x -> x+rc. 
Then we have in (103) in(4+t/*“* 2 ) < £nlx^ I /or x«l and y<<i, so 

tftat corresponding terms can be omitted. The main contriou 
tion to the sum due to the first term ( with unity in the numera- 
tor ) in expression (100) is proportional to 

S(z, 0) = -J yexp(-yz)ln[ l-exp(-y ) ] = lnz/z 2 . (105) 

Thus the number of photons generated in the M-th mode is 

P* = DPmhI 2 * h[ln(M/26)-(-i)*ln(l/26n)]/Otn 2 ) 

N - 1 . _ , , „ 3 n 


(106) 


Since in the considered case (0=2) 6(O=exp(-JTet/L 0 2 )/ir, we get 
‘ 1 'Mowing rate of photon generation for er/Lo »1- 

dP n /dt = 4oS! 0 [1-(-1) h ]/(ji 2 /O. < 107) 


Here Q 0 =^/^o th e mai ' n 

is the dimensionless ampli 
vibrates at the frequency 
for not very large numbers 
arising in approximations 
tion we should limit the t 
sonator t ( due to the diss 
number of photons generate 


eigenfrequency of the resonator . a-e/Lg 
tude of oscillations of the wall ( which 
2 Q a ). Eq. (107) is valid in fact only 
or excized modes M ( due to limitations 
made before). Besides, in real si tua- 
ime t by the relaxation time of the re- 
ipation on the walls). Then the maximum 
d in the M-th mode equals approximately 


n MAX 
r M 


4/(n 2 /0 [2aQ(M)/M + 0(ln(M+i))], 


(108) 


where Q(M) is the quality factor of the resonator' s M-th mode 
This formula is valid provided aQ/M» In (/f+1) . 
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